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Abstract 


This  work  represents  the  joint  work  of  several  groups  including  numerical  calculations  and 
experiments.  Numerical  calculations  were  performed  for  the  propagation  of  a  femtosecond  laser 
beam  through  an  optical  path  and  a  nonlinear  absorber.  The  results  were  compared  with 
experiments  performed  under  similar  conditions.  There  was  excellent  agreement  between 
calculations  and  experiments  at  low  input  energy.  However,  further  additions  must  be  done  to 
the  calculation  of  the  optical  path  for  high  input  energy. 
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1.  Introduction 


Multiphoton  absorption  processes  are  useful  in  several  applications.  These  applications 
include,  blit  are  not  limited  to,  optical  data  storage  [l]-[2],  micro  fabrication  [3]-[7], 
photoconductors  and  pholovoltaics  [8],  markers  for  genomes  and  proteins  [9]-[10],  biological 
and  medical  detectors  [11]-[12],  optica!  limiters  [13]-[16],  biomimetric  electromagnetic  devices 
[17],  nanopattern ing  of  inorganic/organic  materials  [18]-[19]  and  photomedicine/photodynamic 
therapy  [20]. 

The  wide  range  of  applications  of  multiphoton  absorption  processes  makes  it  particularly 
useful  to  have  a  detailed  method  for  numerical  investigations  because  some  of  the  strongest  two- 
photon  absorbers  are  hybrid  chromophores.  These  complex  molecules  exhibit  a  hybrid  photo- 
activated  energy  level  system  in  which  the  two-photon  absorption  (TPA)  level  is  coupled  to  an 
excited  state  absorption  (ESA)  level.  This  hybrid  arrangement  creates  a  complex  dynamical 
system  in  which  the  electron  carrier  concentration  of  every  photo-activated  energy  level  must  be 
taken  into  account  in  order  to  determine  the  actual  optical  properties.  Most  traditional 
calculations  of  the  laser  matter  interaction  make  simplifying  assumptions  about  the  optical  field 
and  the  electron  density  states  of  the  molecular  system.  We  model  propagation  through  two- 
photon  materials  and  describe  the  numerical  analysis  of  the  complex  interaction  of  the  optical 
pulses  with  these  hybrid  systems.  The  numerical  method  calculates  the  spatiotemporal  details  of 
the  electron  population  densities  of  each  photo-activated  energy  level  as  well  as  the  pulse  shape 
in  space  and  time. 

The  main  topic  of  this  effort  was  to  coordinate  a  numerical  simulation  amongst  a 
geographically  and  technically  group  of  experts  (See  list  below).  The  project  involved  a 
simulation  of  an  optical  path  and  propagation  through  a  nonlinear  absorber  with  experiments. 
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Good  agreement  was  found  between  calculations  and  experiments  at  low  input  energy,  but.  at 
very  high  input  energy  it  was  determined  that  additional  term,  such  as  optical  breakdown  in  air 
and  lens,  were  required.  This  report  describes  the  author's  section  of  the  project  involving  the 
propagation  through  the  nonlinear  absorber,  it  necessitated  writing  additional  features  in  the 
propagation  code  including:  (I)  expressing  the  field  in  real  units,  (2)  describing  the  field  in  terms 
of  x,y,z,t  and  r,z,t  and  (3)  inputting  the  field  from  a  digitized  field  file. 

2.  Description  of  the  Test  Case 


The  test  case  involves  simulations  and  comparison  to  experiments.  A  schematic  diagram  of 
the  experimental  set  up  for  the  test  case  is  shown  below. 


1  2  3  4-7  8  9 

Schematic  diagram  of  the  lest  case. 

The  nonlinear  material  is  AF380.  Test  case  parameters:  laser  pulse  width  130  fs  FWHM,  laser 
wavelength  775  nm.  laser  energy  10  microJoule  and  1  mJ,  beam  waist  2.5-3  mm  l/e2  radius, 
divergence  -100  microradians;  sample  thickness)  mm  AF380  film  sandwiched  between  two 
glass  slides,  1  mm  each.  Lens  has  a  focal  length  of  1  m;  lenses  are  either  BK-7  or  SF-1 1  Schott 
glass  with  2.54  cm  diameter.  Data  interchanges  will  be  via  flat  ASCII  files,  SI  units;  in  format 
IpE  14,21;  space  delimited  by  at  least  one  column.  Data  interchange  grid  size:  256x256  (space, 
Cartesian)  x  256  (time)  or  Cylindrical  (256x256  space  x  time).  Field  magnitudes  will  be 
normalized  such  that  peak  amplitude  is  V/m  to  yield  correct  total  energies. 

2.1  Tasks: 
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Table  of  tasks 


Task 

Group 

Activity 

1 

Bill  Dennis  UGA 

model  the  initial  Gaussian  pulse  and  take  care  of  propagation 
through  all  optical  windows/elements. 

2 

Vance  Hedin 

model  all  free-space  propagation  steps  throughout  optical  path. 

3 

Mary  Potasek 

model  propagation  through  all  two-photon  materials  within  setup 

4 

Mark  Walker 

measurements 

5 

ALL 

Final  simulation  results  will  be  passed  back  to  Mark  for  comparison 
to  experimental  data 

Initially  the  data  was  passed  from  one  group  to  another  by  CDs.  Later  a  server  was  set  up  at 
Wright- Patterson  AFB  for  the  data  exchange. 

2.2  Simulation  file  format. 

The  field  file  is  in  the  standard  ASCI!  format  ( I PE24. 1 5E3)  and  order 
(r,t  or  x,y,t).  The  params  file  has  the  following  form: 


Table  for  simulation  file  format. 


Number 

Description 

Units 

128 

Number  of  space  points 

none 

256 

Number  of  time  points. 

none 

1  .OOOOOOOOOOOOOOOE-OO3 

The  width  of  the  spatial  window 

m 

I.474970454238465E-0I2 

The  length  of  the  time  window 

s 

7 . 750000000000000 E-007 

The  center  wavelength 

m 

3.  Numerical  method  used  for  propagation  through  AF380 


The  energy  level  diagram  for  the  AF380  is  given  below. 
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Schematic  energy  level  diagram  of  the  hybrid  multiphoton  absorber  { AF3 80).  The  solid  lines  are  the 
photo  excitations,  the  dashed  lines  show  the  decay  from  the  various  electronic  stales  and  the  dot-dashed 
line  corresponds  to  the  thermal  excitation  directly  into  the  conduction  band  without  any  excited  slate 
absorption.  The  multiple  horizontal  lines  represent  the  vibrational  levels  of  the  various  electronic 
states.The  parameters  are  aipa  =  8.0  e-21  cm4/GW.  alA  =  2.0  e-17  cm2.  <r3  4=2.0  c-17  cm2.  kl(,  =  O.Se-3  ps'1.  ku  =  0.S  c- 
2  ps'1,  kj,j  =  0.29  e-3  ps'1,  lu.j  =  4.0  e-9  ps1,  kj,0  =0.34  e-3  ps'1,  t]  =0.06 

These  parameters  are  used  in  the  rate  equations  to  analyze  all  the  carrier  density  distributions. 
The  propagation  equation  is  coupled  to  these  rate  equations,  which  is  used  to  calculate  the  laser 
pulse  intensity  shape.  Initially  all  of  the  carriers  are  assumed  to  be  in  the  ground  state  (No  ) 
before  the  laser  pulse  enters  the  material.  The  energy  levels  have  the  following  properties.  No  is 
the  ground  state.  N|  is  the  first  excited  state  reached  by  two-photon  absorption.  N;  is  an 
intermediate  state  between  the  first  excited  state  and  the  fluorescing  state  N3.  The  intermediate 
state  is  attributed  to  an  intramolecular  transitional  state  in  which  the  two-photon  excitation 
rearranges  itself  from  the  D-jt  section  of  the  molecule  towards  the  acceptor  region.  Quantum 
mechanical  calculations  of  the  excited  state  electron  density  and  molecular  motion  support  this 
interpretation.  The  upper  level  N4  is  assumed  to  be  the  continuum  (similar  to  a  conduction  band 
in  semiconductors)  that  can  be  populated  from  either  N|  or  N3.  Free  carrier  generation  was 
demonstrated  by  measurement  of  a  nonlinear  photocurrent  using  a  thin  film  of  the  material 
placed  between  transparent  conducting  glasses.  N4  has  a  long  lifetime. 

The  rate  equations  are  given  by 
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where  N,  is  the  electron  number  density  of  the  state  j.  Oj,k  is  the  absorption  cross-section  for 
electron  pumping  from  the  state  j  to  the  state  k,  and  kJ(k  is  the  decay  rate  from  the  state  j  to  the 
state  k.  The  fraction  of  the  two-photon  population  contributing  to  the  photo-induced  current  is 

/  s-1 


given  by  the  Boltzman  distribution  for  electron  rj  = 


1  +  ■ 
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where  T  is  the  temperature,  kn 


is  Boltzman's  constant  and  A E  is  the  energy  gap  between  the  two-photon  state  and  the 
continuum  level. 

The  rate  equations  are  written  in  the  following  form  and  assumed  to  be  in  a  moving 
coordinate  system 

dN  *  —  -  l  -  73-  — 

— —  =  M  N  =[G+——H  +—- — F]N 
dt  tuo0  2ho)0 


where  the  matrix  M  is  given  by 
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The  matrix  M  is  separated  into  three  matrices  for  convenience.  Thus  G  multiplies  the  constant 
decay  terms,  H  multiplies  terms  dependent  on  the  intensity  and  describes  ESA,  and  F  describes 
TPA.  These  matrices  are  given  below. 
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4.  Results 

The  following  section  describes  the  numerical  and  experimental  results.  The  numerical  results 
were  obtained  by  calculating  the  laser  propagation  through  the  lens,  air  and  first  section  of  the 
cuvette  (Prof.  Dennis  group),  followed  by  propagation  through  AF380  (Potasek),  and  then 
propagation  through  the  remainder  of  the  optical  path  (Prof.  Dennis  group).  The  calculations 
were  then  compared  with  experiments  (Mark  Walker  group).  The  following  figures  are  from 
Potasek,  except  where  noted.  The  results  are  divided  into  two  energy  regions:  (I)  low  energy- 
IOuJ  and  (2)  high  energy- 1  OOOuJ.  A  brief  schematic  diagram  of  the  experimental  set  up  is  given 
in  the  figure  below. 


Experimental  Set  up  (Mark  Walker) 


4.1  Low  input  energy-10  uJ 

The  numerical  calculations  and  comparison  with  experiments  are  described  next.  In  this 
case  the  input  energy  is  lOuJ.  The  plots  below  show  the  calculated  contour  plots  for  the  pulse 
propagation  from  the  input  to  the  material  and  the  output  from  the  material.  As  can  be  seen  there 
is  no  pulse  distortion. 


0  rnm  t  mm 


Calculated  contour  plots.  Left  side:  z=0.  Right  side:  z=l  mm. 
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The  figure  below  shows  a  comparison  between  the  numerical  calculation  and  the  experiment.  As 
can  be  from  the  figure,  there  is  good  agreement  between  the  calculation  and  the  experiment. 

Mod/Sim  Comparison 

Time  Integration  of  Input  Pulse  to  Measured  Beam  Profile 
dimensions  are  in  mm 


MLPJ) 

4.2  High  Input  Energy- lOOOuJ 

The  next  section  compares  the  numerical  calculations  with  the  experiments.  The  numerical 
calculations  include  the  calculations  from  both  the  Dennis  group  and  Potasek.  It  should  be  noted 
that  at  high  input  energy  there  is  optical  breakdown  in  the  air  and  cuvette,  which  is  not  taken  into 
account  by  the  Dennis  group  at  this  time.  T  he  plots  below  show  the  calculated  pulse  shapes  and 
contour  plots  as  a  function  of  propagation.  It  can  be  seen  from  the  figures  that  the  pulse  splits  in 
the  temporal  domain.  This  feature  is  not  seen  in  the  low  energy  case. 


0123456789  10 
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measured 


(  From  Mark  Walker. 
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Calculated  Contour  Plots  for  Propagation  through  AF380  for  x  vs  Time 
as  a  Function  of  Distance  (BK-7  lens)  Input  1  mJ 
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Corresponding  Pulse  Shapes  for  x  versus  t  as  a  function  of  propagation  distance. 
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These  temporal  features  are  seen  in  the  comparison  between  calculations  and  experiments  below. 


t 


0  40  60  lJO  !«  200  240 


Instantaneous  beam  profiles  for  x  versus  t.  Left  side:  calculations.  Right  side: 
experiment  (Mark  Walker) 


How'ever.  in  the  xy  pulse  shape  there  is  significant  difference  between  the  calculations  and  the 
experiment.  The  figure  below  shows  the  calculation  on  the  left  and  the  experiment  on  the  right. 
As  can  be  seen  from  the  figure  the  experimental  pulse  shape  is  significantly  larger  than  the 
calculated  one. 


CALCULATION 


9' 
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BK7  IOGOuJ 


Comparison  of  beam  shape  at  high  input  energy.  (Mark  Walker) 
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Additionally,  there  is  a  factor  of  two  difference  between  the  calculated  experimental 
transmission. 

5.  Conclusions 

Comparison  of  numerical  calculations  and  experiments  agreed  well  in  the  case  of  low'  input 
energy.  However,  for  high  input  energy  the  agreement  is  not  as  good.  The  likely  cause  of  the 
difference  between  experiment  and  calculation  is  probably  due  to  neglect  of  ionization  and 
plasma  formation  in  the  air  between  the  lens  and  the  AF380.  Future  calculations  may  need  to 
include  these  effects, 
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7.  Appendix 

Personnel  supported 

PI: 

M.  J.  Potasek 
Collaborators: 

•  Drs.  Sean  Kirkpatrick,  Chris  Brewer  and  Mark  Walker,  WPAFB/ML 

•  Prof.  Dennis,  Univeristy  of  Georgia 

•  Dr.  Erik  Zeek,  Univeristy  of  Georgia 

•  Dr.  Tracey  Bowen,  Kirtland  AFB 

•  Dr.  Vance  Hedin.  Comcast 


Publications: 

No  external  reviewed  publications 
Three  internal  progress  reports 

Interactions/Transitions: 
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a.  Participalion/presentations  at  meeting,  conferences,  etc 

•  Presentation  at  Initial  Meeting.  AFOSR,  Washington.  DC,  Jan  22.  2004 

•  Presentation  at  Review  Meeting,  AT&T,  Fairborn,  Ohio,  Dec  8,  2004 

•  Presentation  at  Review  Meeting,  Kirtland  AFB,  N.M.  Feb.  2,  2006 

•  Numerous  e-mail,  phone  discussions,  and  data  exchange  with  Materials 
Lab,  WPAFB  and  Bill  Dennis  and  Erik  Zeek,  Univ  of  Georgia  and  others 
concerning  the  test  case. 

•  Numerical  code  and  data  exchange  with  Materials  Lab,  WPAFB  and  Bill 
Dennis  and  Erik  Zeek,  Univ  of  Georgia  and  others  concerning  the  test 
case. 

b.  Consultative  and  advisory  functions  to  other  laboratories  and  agencies 

•  Meeting,  discussions,  phone  calls, 

•  Drs.  Sean  Kirkpatrick,  Chris  Brewer  and  Mark  Walker,  WPAFB/ML 

•  Prof.  Dennis,  Univeristy  of  Georgia 

•  Dr.  Erik  Zeek,  Un  iveristy  of  Georgia 

•  Dr.  Tracey  Bowen,  Kirtland  AFB 

•  Dr.  Vance  Hedin,  Comcast 

C.  Transitions 

•  Initial  numerical  code  on  multiphoton  materials  was  transitioned  to  WPAFB. 

New  discoveries,  inventions  or  patent  disclosures: 

None. 
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